The potential for cascading failures in the international trade network.

In our study, we introduce indicators that quantify the influence of each country in complex trade scenarios involving the exchange of raw materials, intermediate goods, and final products across multiple countries. We systematically employ an agent-based model to simulate the propagation of failures from one node to the entire network. This approach allows for the assessment of the impact of each country and the identification of patterns of interaction in the multi-step trade network. Unlike conventional analyses of trade networks, which depict straightforward single-step import/export transactions, our approach captures the intricate realities of processes like raw material procurement, production, and sales in numerous countries from a macroscopic perspective. The findings of our analysis of trade data spanning from 1990 to 2022 reveal several key insights. Firstly, sensitivity to changes in trade volume leading to global failures within interconnected networks has intensified over time. The potential of failure propagation across countries has increased over time, as has the interconnectedness of countries in the global trade landscape. Secondly, despite the increased sensitivity to changes in global trade volume, many countries have become less vulnerable to the influence of others within their multi-step trade networks. This trend aligns with deglobalization, which is evidenced by events such as Brexit and the surge in protectionist measures; these changes indicate a shift in the balance of influence within global trade networks. Thirdly, the results of our analysis of the relationship between load changes and global failures from a regional perspective reveal an intriguing phenomenon: despite limited direct trade connectivity, the interaction between the Latin American and Sub-Saharan African regions is considerable. This suggests the existence of hidden connections between intermediary countries, such that one region's actions can alter the load sensitivity of another, impacting them in unforeseen ways. Furthermore, intra-regional interactions are diminishing in East Asia, while Europe is experiencing a gradual increase in interactions. These trends reflect evolving regional influence, the dynamics of geographic proximity, and the results of economic integration efforts. Additionally, even though the observed period was not long enough to confirm a long-term trend, the previous trend direction was affirmed to persist despite a temporary decrease in trading and reduced sensitivity due to the COVID-19 pandemic. Our study highlights the complexity of global trade dynamics and the need to consider multi-step trade networks and their potential cascading effects when analyzing trade patterns and vulnerabilities.

Porous biodegradable lumbar interbody fusion cage design and fabrication using integrated global-local topology optimization with laser sintering.

Biodegradable cages have received increasing attention for their use in spinal procedures involving interbody fusion to resolve complications associated with the use of nondegradable cages, such as stress shielding and long-term foreign body reaction. However, the relatively weak initial material strength compared to permanent materials and subsequent reduction due to degradation may be problematic. To design a porous biodegradable interbody fusion cage for a preclinical large animal study that can withstand physiological loads while possessing sufficient interconnected porosity for bony bridging and fusion, we developed a multiscale topology optimization technique. Topology optimization at the macroscopic scale provides optimal structural layout that ensures mechanical strength, while optimally designed microstructures, which replace the macroscopic material layout, ensure maximum permeability. Optimally designed cages were fabricated using solid, freeform fabrication of poly(ε-caprolactone) mixed with hydroxyapatite. Compression tests revealed that the yield strength of optimized fusion cages was two times that of typical human lumbar spine loads. Computational analysis further confirmed the mechanical integrity within the human lumbar spine, although the pore structure locally underwent higher stress than yield stress. This optimization technique may be utilized to balance the complex requirements of load-bearing, stress shielding, and interconnected porosity when using biodegradable materials for fusion cages.

A paradigm for the development and evaluation of novel implant topologies for bone fixation: implant design and fabrication.

The future development of bio-integrated devices will improve the functionality of robotic prosthetic limbs. A critical step in the advancement of bio-integrated prostheses will be establishing long-term, secure fixation to the remnant bone. To overcome limitations associated with contemporary bone-anchored prosthetic limbs, we established a paradigm for developing and fabricating novel orthopedic implants undergoing specified loading. A topology optimization scheme was utilized to generate optimal implant macrostructures that minimize deformations near the bone-implant interface. Variations in implant characteristics and interfacial connectivity were investigated to examine how these variables influence the layout of the optimized implant. For enhanced tissue integration, the optimally designed macroscopic geometry of a titanium (Ti)-alloy implant was further modified by introducing optimized microstructures. The complex geometries of selected implants were successfully fabricated using selective laser sintering (SLS) technology. Fabrication accuracy was assessed by comparing volumes and cross-sectional areas of fabricated implants to CAD data. The error of fabricated volume to CAD design volume was less than 8% and differences in cross sectional areas between SEM images of fabricated implants and corresponding cross sections from CAD design were on average less than 9%. We have demonstrated that this computational design method, combined with solid freeform fabrication techniques, provides a versatile way to develop novel orthopedic implants.

Analysis of load sharing on uncovertebral and facet joints at the C5-6 level with implantation of the Bryan, Prestige LP, or ProDisc-C cervical disc prosthesis: an in vivo image-based finite element study.

OBJECT: The goal of this study was to evaluate and compare load sharing of the facet and uncovertebral joints after total cervical disc arthroplasty using 3 different implant designs. METHODS: Three-dimensional voxel finite element models were built for the C5-6 spine unit based on CT images acquired from a candidate patient for cervical disc arthroplasty. Models of facet and uncovertebral joints were added and artificial discs were placed in the intervertebral disc space. Finite element analyses were conducted under normal physiological loads for flexion, extension, and lateral bending to evaluate von Mises stresses and strain energy density (SED) levels at the joints. RESULTS: The Bryan disc imposed the greatest average stress and SED levels at facet and uncovertebral joints with flexion-extension and lateral bending, while the ProDisc-C and Prestige LP discs transferred less load due to their rigid cores. However, all artificial discs showed increased loads at the joints in lateral bending, which may be attributed to direct impinging contact force. CONCLUSIONS: In unconstrained/semiconstrained prostheses with different core rigidity, the shared loads at the joints differ, and greater flexibility may result in greater joint loads. With respect to the 3 artificial discs studied, load sharing of the Bryan disc was highest and was closest to normal load sharing with the facet and uncovertebral joints. The Prestige LP and ProDisc-C carried more load through their rigid core, resulting in decreased load transmission to the facet and uncovertebral joints.

Experimental and computational characterization of designed and fabricated 50:50 PLGA porous scaffolds for human trabecular bone applications.

The present study utilizes image-based computational methods and indirect solid freeform fabrication (SFF) technique to design and fabricate porous scaffolds, and then computationally estimates their elastic modulus and yield stress with experimental validation. 50:50 Poly (lactide-co-glycolide acid) (50:50 PLGA) porous scaffolds were designed using an image-based design technique, fabricated using indirect SFF technique, and characterized using micro-computed tomography (micro-CT) and mechanical testing. Micro-CT data was further used to non-destructively predict the scaffold elastic moduli and yield stress using a voxel-based finite element (FE) method, a technique that could find application in eventual scaffold quality control. Micro-CT data analysis confirmed that the fabricated scaffolds had controlled pore sizes, orthogonally interconnected pores and porosities which were identical to those of the designed files. Mechanical tests revealed that the compressive modulus and yield stresses were in the range of human trabecular bone. The results of FE analysis showed potential stress concentrations inside of the fabricated scaffold due to fabrication defects. Furthermore, the predicted moduli and yield stresses of the FE analysis showed strong correlations with those of the experiments. In the present study, we successfully fabricated scaffolds with designed architectures as well as predicted their mechanical properties in a nondestructive manner.

Topology Optimization of Three Dimensional Tissue Engineering Scaffold Architectures for Prescribed Bulk Modulus and Diffusivity.

Tissue engineering scaffolds play critical roles in skeletal tissue regeneration by supporting physiological loads as well as enhancing cell/tissue migration and formation. These roles can be fulfilled by the functional design of scaffold pore architectures such that the scaffold provides proper mechanical and mass transport environments for new tissue formation. These roles require simultaneous design of mechanical and mass transport properties. In this paper, a numerical homogenization based topology optimization scheme was applied to the design of three dimensional unit microstructures for tissue engineering scaffolds. As measures of mechanical and mass transport environments, target effective bulk modulus and isotropic diffusivity were achieved by optimal design of porous microstructure. Cross property bounds between bulk modulus and diffusivity were adapted to determine feasible design targets for a given porosity. Results demonstrate that designed microstructures could reach cross property bounds for porosity ranging from 30% to 60%.

Stress analysis of the interface between cervical vertebrae end plates and the Bryan, Prestige LP, and ProDisc-C cervical disc prostheses: an in vivo image-based finite element study.

STUDY DESIGN: Segmental motion and bone-implant interface stresses were analyzed at C5-C6 levels with Bryan, Prestige LP, and ProDisc-C cervical disc prostheses using an image-based finite element modeling technique. OBJECTIVE: To predict stress patterns at the interface between prosthesis and lower vertebral end plate to better understand the underlying mechanisms of subsidence and how the load transfer pattern of each disc design affects segmental motion. SUMMARY OF BACKGROUND DATA: Subsidence is one of the most commonly reported device-related complications in intervertebral disc arthroplasty. Although clinical outcomes have been reported regarding many types of cervical prostheses, few reports have analyzed the effects of stress from cervical artificial discs to the vertebral end plate. METHODS: Three-dimensional voxel finite elements were built for C5-C6 spine unit based on computed tomography images acquired from a patient with indication for cervical disc arthroplasty. Models of facet joints and uncovertebral joints were added and artificial disc designs were placed in the intervertebral disc space. Static analyses were conducted under normal physiologic loads in flexion, extension, and lateral bending with precompression. RESULTS: Bryan disc recovered highest range of motion (4.75 degrees ) due to the high elastic nucleus, and therefore imposed the lowest stresses superior to C6. The ProDisc-C and Prestige LP discs caused high stress concentrations around their central fins or teeth, and may initiate bone absorption. Analysis of Prestige LP disc may indicate possible subsidence posteriorly caused by the rear-positioned metal-to-metal joint. CONCLUSION: Rigidity of the cores ("nuclei") in Prestige LP and ProDisc-C prostheses guarantee initial maintenance of disc height, but high contact stress takes place at the bone-end plate interface if they are improperly placed or undersized. Anchorage designs add an additional factor that may increase propensity of subsidence, indicated by the high contact stress occurring at the end plate flanges of Prestige LP, and at midline keel fixation on the end plate of ProDisc-C. Although Bryan disc differs in these 2 concerns, it also creates much larger displacement during motion with more variation in disc height that may theoretically increase the load sharing of facet and/or uncovertebral joints compared to more rigid artificial discs.

Transfer of 90Sr to rice plants after its acute deposition onto flooded paddy soils.

The transfer of 90Sr to rice plants following its acute ground deposition was examined experimentally in a greenhouse. Lysimeters were flooded after being filled with the soil monoliths from 12 paddy fields. A solution of 90Sr was applied to the standing water in the flooded lysimeters at the pre-transplanting stage or booting stage. Applied 90Sr was mixed with the topsoil only after the pre-transplanting application (PTA). The transfer was quantified with the areal transfer factor (TF(a), m2 kg(-1)-dry) defined as the ratio of the plant concentration to the initial ground deposition. In the PTA, the first-year TF(a) values in the 12 soils were in the range of 8.2 x 10(-3) -2.1 x 10(-2) and 1.7 x 10(-4) -3.6 x 10(-4) for the straws and hulled seeds, respectively. The TF(a) values from the booting-stage application (BSA) were higher than those from the PTA by a factor of up to four. The ratios of the seed TF(a) to the straw TF(a) were, on the whole, higher in the BSA. The 90Sr TF(a) in the PTA was negatively correlated with the soil pH and, to a lesser degree, the exchangeable Ca content. In the second year, the TF(a) in the PTA reduced to 53-90% of that in the first year. A more significant reduction, in general, occurred in a sandier soil. Based on the four consecutive years' transfer data, an overall half-time of the 90Sr TF(a) was estimated to be 2.2 years.